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Abstract 

The origin of the Dicke cooperative states, ad hoc assumed for superradiance in the 
system of molecules where no mutual interactions exist but all encountering the same field 
of radiation, was studied by considering two harmonic oscillators driven by a common field 
of radiation. A phasing operator as $/v = D(a)^ P]^D(a) y where Z?(a) is the displacing 
operator and P jv the projection operator for constant energy N for two oscillators, was 
derived. The eigen states of the phasing operator are found to show a finite correlation 
as for the Dicke cooperative states. 


1 Introduction 

The important notion put forward by Dicke [1], that molecules can not be treated as independent 
when the molecules are interacting each with a common field of radiation, introduced the well 
known phenomena of superradiance with the ad hoc assumption of cooperative states. This ad 
hoc assumption of the Dicke’s cooperative states may be a natural consequence of the fact that a 
forced quantum oscillator can be described in terms of the Glauber’s coherent states [2]. However, 
it may be more enlightening to examine in rigorous quantum mechanics how two independent 
molecules (harmonic oscillators) are getting correlated simply by having separate interactions 
with a common field of radiation. 

Furthermore, it may be more interesting if we restrict the interaction between the radiation 
and the molecular system to be “elastic”, that is , no net transfer of energy between the molecular 
system and the radiation field. 

Naively this restriction may correspond to an elastic light scattering from the two molecule 
system and a possibility of phasing or correlation of the two molecules by this continuous scattering 
of light ( more correctly a driving field of radiation in the elastic channel of interaction). 
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2 Two Driven Oscillators 

We consider two harmonic oscillators driven by a classical field : 

H = hu>(af aj + a* 02 ) — xo(a^ + 0.1 + aj + <12 )E(0 


( 1 ) 


where we do not have a direct coupling between the two oscillators. This system seems to be con- 
sidered as trivial because a single harmonic oscillator driven by a classical field is fully understood 

[ 3 ]- 

Since we will be extending the solutions of the single harmonic oscillator driven by a classical 
field to the two oscillator problem of Eq.(l), we want to recollect here some important results of 
a driven harmonic oscillator [3] : 

D(c*)|0) = exp(aa + - a*a)|0) = |a) (2) 

D(a)\n) = ((a + - a*) n /\4d)|a) (3) 

where D(a) is the displacement operator, |a) coherent state of Glauber, 

/ OO 

E(t)expiutdt, exp(aa + - a*a) = exp(-|a| 2 /2) exp(aa + ) exp(-a’a). 

-OO 

We address now to the two oscillator problem of Eq.(l). 

Suppose the two oscillators are prepared in a state | N) of total energy N = ni-fn 2 , then we 
let | N) be driven by a classical field D(o) but we restrict the driven system D(a)|7V) to remain 
at the same energy of |iV). 

The quantum mechanical matrix element corresponding to this process may be written as 


( N\D + (a)P N D(a)\N) 

where Pjv represents the projection operator for states of total energy N : 


N 


Pn = ~ n )l n H n K^ ~ ”1 


n=0 


Alternatively we may define a new operator : 

= D+{a)P N D(a) 

and Eq.(4) can be written as 

(N\* n \N) 

This implies some particular | N) states can become eigen states of the operator $ N . 


(4) 

(5) 

( 6 ) 
(7) 
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3 Evaluation of (iV|$/v|N) 


For two harmonic oscillator states of total energy N we can write 

N 

\N) = ^2°n\N - n)\n) (8) 

n= 0 

From eq.(3) and eq.(8) we obtain 

71=0 

where the suffix (1) and (2) refer to the oscillator index. Equation (5) can be rewritten as 
P N = Y, ((4)"‘ n (4) n Mtf - |0)( 2 )|0)(1){0| (1) {0| (2) (a^~ n /Vn^.^ - n)l) (10) 

71=0 

From eqs.(6), (9) and (10) we obtain 

(JV|4.«|JV) = £ I £ c m (0| (2) (0|(„{«;(«j - 

71=0 771=0 

lyfo>/riy/{N ~ n)\y/{N - ra)!}|a)(i)|<*)(2)| 2 (11) 

Making use of a|or) = a|a) and a n (a + — a*) m = (a + — a*) m a n + m!(a + — a*) m ~ n / (m — n)! 
we can obtain 

N N _ 

(N\$ n \N) = I ^c Tn (0|(2)(0| (1) {v / m!(a^ - a\) m ~ n / Vn\{m - n)\ 

n=0 m=0 

+(4 - a;) m a"/v/^Tv^T}{ y/(N - m)! (a+ - a\) n - m /^(N - n)!(n - m)! 
+(a+ - a *)^-m a ^-n /v / (iV _ m ) !x / (i v _ n)!}|a> (1) |a> (2) | 2 (12) 

Since we are dealing with two identical molecules and the same common field of radiation we may 
set c*i = a 2 = a and make use of the following [3], 

<x|0) = 0, (0|a) = exp(— |a| 2 /2) (13) 


to obtain 

(N\* n \N) 


£ ^^c m c;{exp(-2la| 2 )/n!(Af - n)!}{v / myy / (V - m)! 6 mn 

n=0 m=0 /=0 

+ (-l) N " n >/mI |a| 2 ^- n V(m - n)\y/(N^m)\ 

+ (-\) n ^/(N-m)\ M 2n /(” - m)!Vm! 

+ - ^)! }{y/(N^ly.y/V.6 ln 

+ (-l)"""^ \a\ 2{N ~ n) /{l - n)\\/(N — /)! 

+ (-l)V(^-O! |a| 2n /(” " 

+ (-iy v |a| 2N /v / /! v /(jV - /)! } (14) 
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Since \N) of eq.(8) can be represented by a (N+l)-dimensional state vector, can be defined 
as a (N+l)x(N+l) matrix operator to give 

N N 

(N\i N \N) = Y,'E c ” c '*» o») 

m =0 1=0 

where we find, from eq.(14), the (l,m) matrix element of the operator as follows: 

N 

<b l N = y^(exp — 2|of| 2 /n\(N - n)\){V^.y/(N - m)! S mn 

71 = 0 

+ (-1)<' v -">V^!| 0 | 2 < ,v -”>/( m — n)!v^V — m)! 

+ (-l)V(iV - m)! |a| 2n /(n - m)\Vm\ 

+ (— — ly.Vii s !n 

+ (-1)"-"^/! |a| 2(JV " n) /(f - n)V(JV-/)! 

+ (-i)V(^-O! H 2n /( n ~ 0!v7i 
+ (16) 

We can check for the correct limiting values: 


lirn^(a) 

c*— ►() 


lim (N\$ n \N) 

a— ► 0 


Vim 


N 

E 

71=0 

N 


N 

y Cm (Vm\y/{N - n^./\/n\y/{N - n)lj 6„ 


771=0 


£kl 2 = i 


71=0 


4 <t> N As Phasing Operator 


In order to probe into the physical meaning of the operator we illustrate for the simplest 
nontrivial case of N = ni + = 1 . 

The 4>,v operator is then given in the form 


$N = 1 = 


$oo $01 
$i° $ u 

and the matrix elements are obtained from eq.(16) as follows: 


( 17 ) 


$“ = $}’ = exp(-2|a| 2 )(4 - 8|a| 2 + 8|a| 4 ) 
fc? 1 = $l 0 = exp(-2|a| 2 )(-8|a| 2 + 8H 4 ) 
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Diagonalizing the matrix of eq.(17) gives the eigen values Ai and A 2 as 

A, = $°° + $ 01 , 

A 2 = $°°-$ 01 (18) 

The corresponding eigen states |Aj) and |A 2 ) are obtained as follows: 

I'M = ^/=d°)(l)l 1 )(2) + l 1 )(l)|0)(2)), 

|A») = ^(|0)(i)|l>(j) - |1>(.)|0) W ) (19) 

Correlation or phasing of two oscillators can be measured by a value of < x\ • x 2 > where X\ and 
X 2 are the two harmonic oscillator coordinates of displacement [4]. 

We can see easily 

(0|{2)(l|(l) a: l ' *2|1)(1)|0}(2) = (0|(l)(l|( 2 )*l • Z2|1)(2>|0)(1) = 0 
for the case of a = 0, but for our eigen states of we obtain 

(Aj|xi • x 2 |Ai) = (A 2 |xj • x 2 |A 2 ) 

= {(91(2)0 l(l)**i"^i|9)(i)|l)(j) + (l|<a)(°lu) a i a 2 U)(i)l°}(2)} 

= 4 

We can thus see that the two noninteracting molecules in the common driving field of radiation 
find themselves as correlated. The correlated eigen states of the concerning Hamiltonian of eq.(l) 
can be found as the eigen states of the phasing operator 4>w as introduced in eq.(6). 

The existence of the correlated eigen states of the two oscillator Hamiltonian thus justifies 
the ad hoc assumption of the Dicke cooperative states for independent molecules all in the same 
common field of radiation. 
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